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The study of the Higgs couplings following its discovery is the priority of future LHC runs. A hint 
of anomalous nature will be exhibited via its coupling to the Standard Model(SM) particles and open 
up new domain of phenomenological study of physics beyond the Standard Model. The enhanced 
statistics from next LHC runs will enable entry into the precision era to study the properties of 
Higgs with greater details. In this paper we present how one can extract Higgs couplings in future 
LHC runs at 14 TeV via H —>• ZZ* —> 4 1, using observables constructed from angular distributions 
for the Standard Model Higgs and Higgs with mixed CP configuration. We show how angular 
asymmetries can be used to measure the ratios of the couplings and the relative phases at LHC. We 
benchmark our analysis finding out the angular asymmetries and the best fit values of the ratios of 
the couplings for SM Higgs, CP-odd admixture, CP-even higher derivative contribution and when 
CP-even higher derivative contribution and CP-odd admixture are both present. In the Standard 
Model, HZZ couplings have no momentum dependence. It is thus essential to demonstrate the 
momentum independence of the couplings to establish the couplings are SM like in nature. In this 
work we show how one can test the momentum independence of the Standard Model like coupling 
using angular asymmetries. We develop the necessary tools and demonstrate how to study the 
momentum dependence can be studied at future LHC runs. 


I. INTRODUCTION 

The ATLAS [1, 2] and CMS [3-5] collaboration at LHC 
have both discovered a new resonance of mass around 
125 GeV that is found to be largely consistent with 
the observation of a Higgs boson. Several studies [6- 
59] done both before and after the discovery of the Higgs 
boson have examined how to determine the spin, parity 
and coupling of the Higgs boson. In gauge sector, decay 
modes such as H —> 77, H —> ZZ and H —> WW etc., 
where one (or both) of the Z’s and W’s are off-shell, are 
used to study the spin, parity and coupling of the Higgs 
boson. Observation of the decay mode H —> 77 estab¬ 
lishes that the discovered resonance is necessarily a boson 
and the Landau-Yang theorem [65, 66] excludes the possi¬ 
bility of it having spin J = 1. Furthermore if Higgs boson 
is a eigenstate of charge conjugation, charge conjugation 
invariance along with observation of H —> 77 also en¬ 
force [10] that Higgs is a charge conjugation C = + state. 
Recent measurements [70-73] have shown that the reso¬ 
nance favors Spin 0 over spin 2. Moreover Ref. [55] rules 
out pure pseudoscalar hypothesis i.e. J p = 0“ at a 99.98 
% CL. However the discovered Higgs can still have small 
CP-odd admixture or higher derivative CP-even contri¬ 
bution to its coupling. Angular distributions and angular 
asymmetries of Higgs decay are essential to investigate 
whether the discovered resonance is a CP eigenstate or 
a resonance with mixed CP configuration. As these an¬ 
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gular asymmetries are functions of the Higgs couplings, 
studying them will allow us to probe the nature of the 
Higgs coupling directly. 

In this paper we will restrict ourselves to experimen¬ 
tally clean golden channel H —*• ZZ* 
where I\, I 2 are leptons e or /i. We consider Spin-0 Higgs 
boson H with even parity but, include the possibility of a 
small CP-odd admixture and higher derivative CP-even 
contribution. We first calculate the differential decay rate 
for H — > ZZ* — > (^+)(4-£+) process in terms of in¬ 
variant mass of the dileptons coming from the off shell Z 
boson and angular distributions of the final state leptons. 
From these distributions, we construct angular asymme¬ 
tries (observables) and utilize them to probe anomaly in 
HZZ couplings. Similar asymmetries have also been dis¬ 
cussed in Ref . [60-62] . As these observables are functions 
of HZZ vertex factors, the values of the different observ¬ 
ables differ for the various cases such as SM Higgs, CP- 
odd admixture, CP-even higher derivative contribution. 
Ref. [53, 55] have discussed how ratios of couplings can 
be measured at 8 TeV LHC run. Ref. [53] shows how us¬ 
ing uniangular distribution as input in likelihood analysis 
one can discriminate different spin possibilities. In this 
paper a study showing how using simple angular asym¬ 
metries one can study the CP property of H in future 
LHC runs. 

We benchmark these observables for SM Higgs and 
Higgs with a CP-odd admixture at 14 TeV 300 fb” * 1 LHC. 
We determine the ratios of the coupling constants and use 
them to discriminate possible CP-odd admixture from 
SM Higgs. We then perform the same analysis at 3000 
fb' 1 to distinguish SM, case with CP-odd admixture and 
CP-even higher derivative contribution. We also consider 
the scenario when higher derivative CP-even and and CP- 
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odd admixture are present in HZZ couplings. We denote 
this scenario as ‘CP-even-odd’ case. In our analysis we 
have also included a complex phase for CP-odd admix¬ 
ture Higgs and show how to determine the phase using 
angular asymmetries. Furthermore we use these angular 
asymmetries and perform chi-square analysis to probe 
both CP-even and CP-odd anomalous contributions in 
the Higgs couplings. 

Exotic models of Higgs can have momentum depen¬ 
dence in its couplings. It is thus essential to study the 
momentum dependence of the Higgs couplings to estab¬ 
lish its SM nature. Angular asymmetries will provide 
necessary tools to investigate the momentum dependence 
of the Higgs couplings. In our work we develop necessary 
techniques and utilize them to probe the momentum de¬ 
pendence of HZZ couplings. It is important to note this 
it will definitely require higher statistics. However the 
enhanced statistics at 14 TeV, LHC will enable us to 
study the momentum dependence of Higgs couplings in 
different momentum regions. Since the mass of H does 
not allow both the Z bosons to be on-shell, the invari¬ 
ant mass distribution of the dileptons from the off-shell 
Z boson(Z 2 ) will offer us a test for the momentum de¬ 
pendence of the HZZ couplings. The HZZ couplings in 
the most general case could be function of the invariant 
mass of the off-shell Z. The constancy of the ratios of 
the Higgs couplings can be measured by finding out the 
values of the ratios of couplings in different momentum 
regions for the invariant mass of the off-shell Z. In our 
analysis we have shown for SM how one can test the mo¬ 
mentum dependence of the Higgs couplings at 14 TeV 
300 fb -1 LHC run. 

In LHC we measure a ■ BR and that does not allow 
us to measure the decay width of the Higgs, as a re¬ 
sult we can only measure the ratios of Higgs couplings. 
However ILC[63] will be able to measure the inclusive 
cross section(crzjy) using recoil technique for the process 
e + e _ —> ZH. Hence inclusive cross section provide a di¬ 
rect measurement of HZZ couplings at ILC by measuring 
partial width T(H ZZ). After energy upgrade LHC 
will run at 33 TeV and enhanced statistics will enable 
us to measure Higgs couplings more precisely than that 
possible at 14 TeV. We discuss how much this energy up¬ 
grade will improve the measurement of HZZ couplings 
and discuss what sensitivity can be achieved compared 
to that of an ILC[64] measurement. 

The paper is organized as follows: In Section II we for¬ 
malize the necessary tools for our analysis. Section III 
is divided into three subsection. First we benchmark 
our analysis for 300 fb -1 in Subsection III A and then 
we benchmark our analysis for 3000 fb -1 for SM Higgs, 
Higgs with CP-odd admixture, CP-even higher deriva¬ 
tive contribution and CP-even-odd scenario in Subsec¬ 
tion IIIB. In Subsection IIIC we develop the neces¬ 
sary technique to measure the momentum dependence 
of Higgs couplings. A qualitative comparison between 
future runs of LHC and ILC precision measurement has 
been made in Section IV. Finally we conclude in Sec¬ 


tion V. 


II. THE FORMALISM 


In this section we first write down the HZZ vertex, 
helicity amplitudes in transversity basis and finally derive 
the expression for angular distribution of H —> ZZ* —> M 
process assuming H to be a spin 0 particle. In SM the 
process H —► ZZ is characterised by the Lagrangian 


£-hzz 


gMz Z p Z p H 
2 cos 9w 


(1) 


where 6\y is the Weinberg angle and g is the electroweak 
coupling constant. However there may exist anomalous 
couplings of H to Z boson. These couplings can in gen¬ 
eral be CP-even or CP-odd and can be generated from 
the effective Lagrangians 

£ e ~ —Z^Z^H ( 2 ) 

and 

Co^—Z^Z^H (3) 


respectively, where Z pv and Z^ are defined as Z= 
d l ,Z u - d v Z p and Z pv = e pupa Z p ° respectively. 

Following these Lagrangians one can write down the 
most general HZZ vertex as follows 

V^ = ^r(agr + b (<#<£ - qi -q 2 gT) 

cost 1\v \ 

+ ic t pvpa qip q 2 a' S J , (4) 


where a, b , c are momentum dependent vertex factors 
and < 7 i, q 2 and P are the four momenta of Z(Z i), Z*(Z 2 ) 
and H respectively. Off-shcllness of the Z is denoted by 
the superscript In Standard Model at tree level the 
values of the vertex factors are a = 1 and b = c = 0 
and they are constant. However a non zero b and c can 
arise from higher order correction. If the vertex factors 
a , &, c show any deviations from SM values or exhibit 
a momentum dependence it would provide a hint about 
the non standard nature of the HZZ couplings. The CP- 
ood admixture is charactersided by the non zero value of 
c of the form c e lS , where <5 is the CP violating phase 
associated with c. 

The decay under consideration can be characterised by 
three helicity amplitudes Al, A\\ and A± defined as (see 
appendix) : 

Al = qi ■ q 2 (a - b q 1 .q 2 ) + M% X 2 b, (5) 
A\\ = \J%qlql ( a-b qi .q 2 ) , ( 6 ) 

A± = y/2qlq*XM H c, (7) 
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where y/qf and yfq f are the invariant masses of the Z\ 
and Zi respectively, with 


\J 1,92) 

2 M h 


( 8 ) 


and 


X(x, y,z)=x 2 +y 2 + z 2 — 2 xy— 2 xz — 2 yz. (9) 


FIG. 1. The definition of the polar angles ( 6 1 and 62 ) and 
the azimuthal angle (< t >) in the decay of Higgs (H) to a pair 
of Z boson, followed by a decay of both the Z into four 
charged leptons: H —> Z\ + Z 2 —> (if + if) + (if + if)-, 
where l\,l 2 £ {e, y} and three momentum k\ = — and 
£3 = — £ 4 .The lepton pair if going back to back in the rest 
frame of Z \, however lepton pair from Z 2 going back to back 
in their center-of-momentum (C.O.M) frame. 


The helicity amplitudes have definite parity properties as 
Al, A\\ are CP-even and A± is CP-odd. 

Having defined the helicity amplitudes in transversity 
basis the full angular distribution for H —> Z\ + Z 2 —> 
(£f+£+) + (£f+i+), can be written as [53] 


7 -ffTf -5—7-7 ) 3T = !+ '-7- cos - P 2 (cos6»i))( 1 -P 2 (cos0 2 )) + -Im(.F||.FjO sin 2(f) 

If UQ 2 a costal a cos 62 dcp 4 2 " 

X (l - P 2 (cos0i))(l - P 2 (cos0 2 )) + ^(1 - 3 IJ'lI 2 ) (P 2 (cos0i) + P 2 (cos0 2 )) + ^(1 + 3 |Phi 2 ) P 2 (cos0i)P 2 (cos0 2 ) 

—^Re(J r LJ r |*) cos <j> + Im(PiPp) sinpj sin 20i sin 20 2 + 77^2_Re(J r LJ r ][)(cos0i — cos0 2 ) cos<^>sin0i sin0 2 

9 

— cos 0i (2 + P 2 (cos0 2 ))) — ry^-^=Im(J r iJ r |*)( cos0i — cos0 2 ) sin0sin0i sin0 2 

( 10 ) 


COS (j>- 


+ 8y/2 

+ ^Re(P||.Fl)(cos0 2 (2 + P 2 (cos0i)) — cos0i(2 + P 2 (cos0 2 ))) - ty^=Im(7' i 7'p(coi 
- ^r ? 2 ^(1 - |Pz, | 2 ) cos 0 i cos 0 2 + V 2 (Re(Pz,Pp cos </> + Im(Pi,Pp) sin (fj sin 0 i sin 0 2 


where Ph, P| i, Pj _ are the helicity fractions defined in 
the appendix along with Tf and rj. The angle 0i(0 2 ) is 
the angle between three momenta of If (If) in Z\(Zf) 
rest frame and the direction of three momenta of Z\(Zf) 
in H rest frame. The angle <f> is defined as the angle be¬ 
tween the normals to the planes defined by Z\ —► if if 
and i? 2 —► if if in H rest frame as shown in Fig. 1. It 
should be noted that Eq.(10) is exact and no assump¬ 
tions has been made apart from assuming that the leptons 


are massless. 

Integrating Eq.(10) with respect to any two out of the 
three angles 0 i, 0 2 , <f> one finds three uniangular distri¬ 
butions for the process H ZZ —> (ifif)(ifif) as: 

77— , 2 ^ r a = ^ +r 2 P 2 (cOS0i) -PCOS01, (11) 
Tf dq 2 d cos 0i 2 

77 - , 2 ^ r a = \ +T 2 P 2 (cos 0 2 ) + Tx cos 0 2 , ( 12 ) 

If aqZ a cos 62 2 
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— , .) , . = 1 + C /2 cos 2 <p + V 2 sin 2</> 

( 2^2 CL(p 

+ U\ cos 4> + Vi sin (13) 

The observables Ti, T 2 , f7i, (72, Vi and V 2 are defined 
as: 


Ti = ^Re(T||TI), 

(14) 

T 2 = i(1^3|T L | 2 ), 

(15) 

Q7T 2 

(16) 

U 2 = J(|T I || 2 -|T]_| 2 ), 

(17) 

Qtt 2 

y ‘ = -32V2” 

(18) 

U 2 = ilm(T||Pl), 

(19) 


and are functions of q§. Moreover Pi (cos $i j2 ) P 2 (cos $i :2 ) 
are first and second degree Legendre Polynomials respec¬ 
tively. It should be noted that the observables Ti, P 2 , 
Ui, U 2 , V\ and V 2 are coefficients of orthogonal func¬ 


tions P 2 (cos # 1 , 2 ), Pi(cos0i j2 ), cos2 (j>, cos</>, sin20, sine/) 
respectively and can be extracted individually. These 
observables are functions of the helicity fractions Tl, F\\ 
and Tj_, hence they are functions of vertex factors a, b 
and c. Measurement of these observables will enable us 
to probe the vertex factors a, b and c. Moreover in SM 
T 2 , U 2 and U\ all are non zero as Tl and P) are non 
zero. In SM at tree level c = 0 which enforces T_\_ = 0. 
As Ti, V 2 and Vi all are functions of Pj_ = 0, in SM they 
are all zero. A CP-odd admixture is characterised by 
non zero value of c and hence non vanishing values of T \, 
V 2 and Vi. Measurements of Ti, V 2 and Vi allow us to 
probe CP violating phase in HZZ couplings. If there ex¬ 
ist any CP-even higher derivative contribution in HZZ 
couplings, the vertex factor b becomes non zero, hence 
the observables T 2 , (7 2 and U\ will have different values 
than that of SM. In SM at tree level b = 0 but at one loop 
level the value of b will be non zero. Measurement of b 
will allow us to probe triple-Higgs vertex which arises at 
one loop level and provide the first verification of Higgs 
self coupling. 

This observables can be extracted using following an¬ 
gular asymmetries 


Ti 




1 d 2 r \ 

Tf dq% d cos ) 





U, = 


r+D 


'+h 


dcos0i j2 


1 


d 2 r 




(* + 7T > 


'+'9 


d(j) 


Pf dcos0i j2 
2tt d 2 T 


If dq 2 d(f> 


- ;.(// /:•/• r/; 


d(f> 


d 2 T 
dq% d<j> ’ 



1 d 2 r \ 

Tf d cos d 2 ) ’ 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 


Three uniangular distributions in Eq.(ll), Eq.(12) and 
Eq.(13) give 6 observables Ti, T 2 , (7i, V\ and V 2 in terms 
of angular asymmetries shown above. At 14 TeV, high 
luminosity future LHC runs will provide enhanced statis¬ 
tics and the uniangular distributions will not only allow 
us to probe Higgs couplings, but also provide us with ar¬ 
senal to probe the momentum dependence of the Higgs 
couplings which is essential for precision measurement. 
In the next section we will discuss these aspects in de¬ 
tail. 


III. NUMERICAL STUDY 


In this section we demonstrate how angular analy¬ 
sis can be used to benchmark different CP scenarios of 
H. We have generated events with MADEVENT5 [75] 
event generator interfaced with PYTHIA6.4 [76] and 
Delphes 3 [77]. The vertex, Eq. (4) is parametrized by 
UFO format of MadGraph5 using HiggsCharacterisation 
model [78]. The events are generated by pp collisions 
via gg —> H and gg H + ljet, for center of mass en- 
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ergy yfs = 14 TeV, using parton distribution functions 
CTEQ6L1 [74], 

For matching purpose we have used MLM prescription 
and events are finally passed through fast detector sim¬ 
ulator package Delphes 3. We are only concerned about 
the decay process of the H and have assumed only the 
SM production process of H while generating events. 

We follow the cut based analysis of Ref. [73]. The iden¬ 
tical pairs of final state leptons 2 e + 2 e _ and 2 p + 2p— 
events have also been taken into account for our analyses. 
The on-shell Z\ boson is identified by the invariant mass 
of the opposite charge same flavor lepton pairs closest to 
Mz- Moreover, since the Higgs mass Mh does not al¬ 
low both the Z bosons to be on-shell, this in turn breaks 
the need for Fermi antisymmetrisation of the identical 
fermions when the final states are 2 e + 2 e - and 2 /i + 2 /x—. 
The branching ratios and cross sections have been taken 
from Higgs working Group webpage [79]. 

Following the analysis presented in Ref. [73] data are 
selected using single-lepton or di-lepton triggers. For 
the single-muon trigger the transverse momentum, Pt , 
threshold is 25 GeV, while for single-electron trigger the 
transverse energy, Et, threshold is 25 GeV. Di-muon 
triggers are selected using two ways. For asymmetric di¬ 
muon the trigger thresholds are either pri = 18 GeV and 
PT 2 = 8 GeV. Threshold for symmetric di-muon triggers 
are pr = 13 GeV for both the muons. For the di-electron 
trigger the thresholds are Et = 12 GeV for both elec¬ 
trons. There are two electron-muon triggers used with 
12 or 24 GeV Et electron thresholds, differed by the 
electron identification requirements, and muon threshold 
p T = 8 GeV. 

Each electron (muon) must satisfy Et > 7 GeV (pr > 
6 GeV) and be measured in the pseudo-rapidity range 
1 77 1 < 2.47 (| p | < 2.7). We have selected the leptons in 
two sequential pt ordered way. 

i) Case-I : pt of at least two leptons in a quadruplet 
must satisfy pr > 20 GeV, 

ii) Case-II : pt of at least three leptons in a quadruplet 
must satisfy pt > 20 GeV. 


The leptons are required to be separated from each 
other by A R > 0.1 if they are of the same flavour and 
A R > 0.2 otherwise. Each event is required to have the 
triggering lepton(s) correctly matched to one or two of 
the selected leptons. 

Furthermore we also impose the invariant mass cuts on 
the mz 2 {\/~q f) and mu described in Table I. 

is the invariant mass of the pair of opposite sign 
same flavor leptons closest to mz while mz 2 is the other 
combination. The two columns of Table I demonstrate 
the effect of pt ordering in event selection. 

Now integrating Eq. (11),(12), (13) over qf we get three 
integrated distributions as follows 

1 dr _ 1 

r d cos 61 2 


Cuts 

Case-I 

Case-II 

Selection cuts 

494 

2253 

50 GeV < mi 2 < 106 GeV 

487 

2204 

12 GeV < 77134 <115 GeV 

447 

2071 

115 GeV < mu < 130 GeV 

443 

2050 


TABLE I. Effects of the sequential cuts on the simulated Sig¬ 
nal for two different pr ordering of Case-I(first column) and 
Case-II( second column). The sequential pr ordering of Case- 
I is for 300 fb -1 , however we have used sequential pr ordering 
of Case-II for 3000 fb -1 . The k -factor for signal is 2.5. 


= 7 : +71 cos 6> 2 +7^P 2 (cos6* 2 ), (27) 

1 CL COS U2 ^ 

4 = 77- + U\ cos (f> + U 2 cos 2 <j> 

1 dtp 27t 

+ Vi sin <j>-1 - V2 sin 2 </>, (28) 


where 71, 72, U\, U 2 , Vi and V 2 are observables inte¬ 
grated over 77134 ( 52 ) and 77112 . 

The normalized distributions, YJcoTe^ vs cos ^i > 
T d cos e 2 vs cos ^2 and vs <j) for SM are shown in 
Fig. 2, Fig. 3 and Fig. 4 respectively for simulated data. 
It should be noted that the angular coverage for cos 61 or 
cos 62 covers the full range from —1 to +1 and coverage 
for (f> from 0 to 27 t are still retained even after using ac¬ 
tual detector scenarios. The cut flow analysis of Case-I 
is followed for the analysis of SM Higgs and Higgs with 
CP-odd admixture at 300 fb -1 . At 3000 fb -1 since the 
statistics is higher, we will use stronger cut based analysis 
i.e. sequential cut flow analysis of Case-II for benchmark¬ 
ing SM Higgs and Higgs with different CP configuration. 
Moreover it should be noted that we have used the same 
cut based analysis for CP-odd admixture, CP-even higher 
derivative contribution and CP-even-odd scenario. The 
cross section for each benchmark scenarios are within the 
current experimental allowed region. 


-1 -0.8 -0.6 .0.4 -0.2 


simulated result - 


0.2 0.4 0.6 0.8 1 


FIG. 2. The normalized distribution ^ . dr „ vs cos#i for SM 

r a cos 6 1 A 

Higgs events. 


71 cos#i + 72 P 2 (cos#i), (26) 
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simulated result - 


-1 - 0.8 - 0.6 - 0.4 - 0.2 0 0.2 0.4 0.6 0.8 1 


FIG. 3. The normalized distribution -k , dl a vs cos 62 for SM 

I a cos 02 

Higgs events. 


% 0.16 - 


simulated result - 


FIG. 4. The normalized distribution vs <j> for SM Higgs 
events. 


The simulated data are binned in cos 6 b, cos 62 and </> 
and fitted using Eq.(26), Eq.(27) and Eq.(28) to obtain 
the angular asymmetries t\, 62 , mi, M 2 , v\, V 2 and their 
errors which correspond to the angular asymmetries 71, 
71, Hi, U 2 , Vi and V 2 respectively. Once the values of 
the integrated observables t\, t. 2 , Mi, M 2 , v\, V 2 and their 
respective errors are found, the x 2 formula: 

2 _ (72 t 2 ) 2 (cos S 71 H) 2 (W 1 -M 1) 2 

X ~ (At 2 ) 2 (Afr) 2 + (Ami) 2 

{U 2 —U 2) 2 (sin 5 V 2 — m 2 ) 2 (sin (5 Vi — Mi ) 2 
(AM 2) 2 (AM 2) 2 (Ami) 2 

(29) 

will find the b/a, c/a and the phase S. The errors in b/a, 
c/a and phase 6 can also be calculated using the error 
matrix 


d 2 X 2 \ 
daidaij - 


( 30 ) 


where ai,atj = b/a, c/a, 5. To find the best fit values we 
have used Mathematica 9 [69]. 


A. Study of angular asymmetries of Higgs at 14 
TeV and 300 fb 1 

We start benchmarking angular observables for SM 
Higgs and Higgs with CP-odd admixture in this section. 
The fit values for observables are tabulated along with 
the best fit values of b/a, c/a and phase S of the CP- 
odd coupling. We also obtain 1 <t and 2er contours for b/a 
vs c/a and S vs c /a for the case of CP-odd admixture. 
This will provide the precision at which one can rule out 
anomalous contributions in HZZ couplings, establishing 
the SM nature of H at 14 TeV 300 fb -1 LHC run. 


1. SM Higgs 

SM Higgs events are generated with a = 1, b = 0 , 
c = 0. The fit values of the observables for the SM Higgs 
are tabulated in Table II. 

TABLE II. The values of the observables for the SM Higgs 
with respective errors. 


Observables 

Values with errors 

£2 

-0.21 ± 0.09 

ti 

(-1.6 ±7.16) x 10' 2 

U2 

0.32 ± 0.40 

Ml 

(0.93 ±4.36) x 10" 1 

V2 

(-0.72 ±4.03) x 10 _1 

Ml 

(0.19 ±3.70) x 10" 1 


The values of the observables t 2 and M 2 are large com¬ 
pared to other observables as discussed in the previous 
section, playing important role in the y 2 expression in 
Eq.(29). The observables t\, V 2 and mi provide informa¬ 
tion about phase for anomalous couplings b and c. The 
best fit values of b/a and c/a with their respective errors 
for the SM Higgs are given as follows: 

b/a = (0.50 ± 0.96) x 10 ” 4 GeV “ 2 (31) 

c/a = (0.68 ± 2.27) x 10 ~ 4 GeV “ 2 (32) 

The best fit values with la and 2<r contours for b/a vs 
c/a are shown in Fig .5. 


2. Higgs with CP-odd admixture 

The CP-odd admixture is charaterised by a non zero 
value of c in Eq.(4). For CP-odd admixture case, Higgs 
events are generated using a = 0.7, 6 = 0 and c = (2.2 + 
2.2 i) x 10 -4 . The values of the observables are given in 
Table III. 



















































7 



FIG. 5. c/a vs b/a la (green) and 2a (yellow) contours for 
the SM Higgs at 300 fb -1 . The best fit values (b/a, c/a) is 
shown by the block dot. The V corresponds to b = c = 0. 


TABLE III. The values of the observables for CP-odd admix¬ 
ture Higgs with respective errors at 300 fb -1 . 


Observables 

Values with errors 

t 2 

-0.06 + 0.10 

ti 

-0.11 + 0.08 

U2 

-0.08 ± 0.40 

Ui 

(-0.80 ± 4.04) x 10" 1 

Vi 

(0.99 + 4.20) x 10 _1 

Vi 

(0.44 + 4.21) x 10 1 


The value of f 2 has now become smaller compared to 
the SM case as shown in Table II. Most importantly the 
non zero value of + arises due to the complex CP-odd 
anomalous coupling c. This will play a significant role 
along with f 2 and m 2 in probing anomalous CP-odd ad¬ 
mixture of HZZ couplings. The best fit values for b/a, 
c/a and the phase 8 for CP-odd admixture are : 

b/a = (1.50 ± 1.09) x 10' 4 GeV“ 2 (33) 

c/a = (5.48 ± 1.12) x 10” 4 GeV“ 2 (34) 

d = (0.29 ± 2.14) in radian. (35) 

Note that the error in 8 is still very large at this lumi¬ 
nosity. 

The best fit values with lcr and 2cr contours for c/a vs 
b/a and 8 vs c/a are shown in Fig. 6 and Fig. 7 respec¬ 
tively. 


FIG. 6. c/a vs b/a la (green) and 2cr (yellow) contours for 
CP-odd admixture Higgs at 300 fb -1 . The best fit value of 
(b/a, c/a) is shown by the block dot. The values with which 
data are generated (b/a = 0, c/a = 4.44 x 10 -4 ) is shown by 
the ‘> 1 =’. The cross-hair corresponds to b = c = 0. 



c/a (10“ 4 GeV -2 ) 

FIG. 7. 5 vs c/a la (green) and 2cr (yellow) contours for CP- 
odd admixture Higgs at 300 fb -1 . The best fit values (c/a, 8) 
is shown by the block dot. The values with which data are 
generated is shown by the V. 

3000 fb -1 also, we have followed the same cut based anal¬ 
ysis that we have discussed earlier apart from a strong 
sequential px ordering i.e. px of at least three leptons in 
a quadruplet must satisfy px > 20 GeV. 

At 3000 fb -1 we revisit the benchmark cases of SM 
and CP-odd admixture along with two new analysis of 
CP-even higher derivative contribution and CP-even-odd 
scenario. 


B. Study of angular asymmetries of Higgs at 14 
TeV 3000 fb 1 

High Luminosity LHC (HL-LHC) i.e 14 TeV 3000 fb -1 
run before the energy upgrade will allow us to test CP 
structure of HZZ couplings even more precisely. For 


1. SM Higgs and CP-odd admixture Higgs 

First we investigate CP-odd Higgs and SM Higgs and 
find out the values of angular observables along with their 
respective errors. For CP-odd admixture we have again 
taken a = 0.7, b = 0, c = (2.2+2.2*) x 10 -4 and SM Higgs 


















a=l,6 = 0,c=0. The fit values of the observables £ 2 , 
£ 1 , U 2 , ui, V 2 , v\ for the SM and CP-odd admixture Higgs 
are tabulated in Table IV and Table V respectively. The 

TABLE IV. The values of the observables for the SM Higgs 
with respective errors at 3000 fb” 1 . 


Observables 

Values with errors 

^2 

-0.20 ± 0.04 

ti 

(0.28 ±0.35) x 10” 1 

U2 

0.21 ±0.19 

Ml 

(0.46 ±2.06) x 10” 1 

V2 

(-0.07 ± 1.96) x 10” 1 

Ml 

(-0.18 ± 1.84) x 10” 1 


errors have significantly reduced for all the observables 
and the fit values for ratios of couplings for the SM Higgs 
b/a, c/a are given 

b/a = (0.43 ± 0.55) x 10” 4 GeV” 2 (36) 

c/a = (1.08 ± 1.17) x 10” 4 GeV” 2 (37) 

The lcr and 2<r contours for b/a vs c/a are shown in 
Fig. 8 



FIG. 8. c/a vs b/a la (green) and 2cr (yellow) contours for 
the SM Higgs at 3000 fb -1 . The best fit values (b/a, c/a) is 
shown by the block dot. The V corresponds to b = c = 0. 

At 3000 fb -1 from Table V one can see that the errors 
in t\ and £2 are much reduced, making them very good 
observables for probing CP-odd admixture. The best fit 
values for b/a, c/a and phase 5 for CP-odd admixture 
are given as 

b/a = (0.40 ± 0.54) x 10” 4 GeV” 2 (38) 

c/a = (3.99 ± 0.64) x 10” 4 GeV -2 (39) 

5 = 0.45 ± 1.08 in radian. (40) 

It should be noted that the error in 5 has become lower 
due to the fact that the error in t \, which constraints the 
phase S, is much reduced. 


TABLE V. The values of the observables for CP-odd admix¬ 
ture Higgs with respective errors at 3000 fb -1 


Observables 

Values with errors 

t2 

-0.11 ± 0.04 

ti 

-0.06 ± 0.03 

U2 

0.02 ±0.18 

Ml 

(-0.10 ±0.56) x 10” 1 

V 2 

(0.72 ± 1.84) x 10” 1 

Ml 

(0.67 ± 1.83) x 10” 1 


The lcr and 2er contours for b/a vs c/a and 6 vs c/a 
are shown in Fig. 9 and Fig. 10 respectively. 



FIG. 9. c/a vs b/a la (green) and 2cr (yellow) contours for 
CP-odd admixture Higgs at 3000 fb” 1 . The best fit value of 
(b/a, c/a) is shown by the block dot. The values with which 
data are generated (b/a = 0, c/a = 4.44 x 10” 4 ) is shown by 
the V. The cross-hair corresponds to b = c = 0. 

So far, we have discussed how using angular asymme¬ 
tries one can probe HZZ couplings of the SM Higgs and 
Higgs with mixed CP scenarios at 14 TeV for two differ¬ 
ent luminosity 300 fb” 1 and 3000 fb” 1 . The values of the 
observables vary depending on the values of a, b and c. 
The observables Ti, Vi and V 2 are sensitive to CP-odd 
admixture and can be a good candidate to probe CP-odd 
admixture. Finally the best fit values for the ratios of the 
couplings, b/a and c/a are calculated using Eq.(29). 

2. Higgs with CP-even higher derivative contribution and 
CP-even-odd scenario 

In this subsection we benchmark the angular asymme¬ 
tries for the cases: 

1) Higgs with CP-even higher derivative contribution. 

2) Higgs with both CP-odd and CP-even higher deriva¬ 
tive contribution i.e. CP-even-odd scenario. 

For CP-even higher derivative contribution we have taken 
a = 0.80, b = 10” 4 , c = 0. However for the CP-even- 
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FIG. 10. 5 vs c/a la (green) and 2cr (yellow) contours for 

CP-odd admixture Higgs at 3000 fb -1 . The best fit values 
(c/a, 8) is shown by the block dot. The cross-hair corresponds 
to 5 = c = 0. The value with which data are generated is 
shown by the V. 

odd scenario we have taken a = 0.75, b = 8 x 10 -5 , 
c = (1 ± li) x 1CT 4 . 

The values of the observables for CP-even higher 
derivative contribution are tabulated in Table VI. The 

TABLE VI. The values of the observables for Higgs with CP- 
even higher derivative contribution and their respective errors 
at 3000 fb" 1 . 


Observables 

Values with errors 

£2 

-0.12 ±0.04 

ti 

(0.31 ±3.40) x 10" 2 

«2 

0.16 ±0.18 

Ul 

(-0.2 ± 1.90) x 10" 1 

Vi 

(-0.09 ± 1.82) x 10" 1 

Vi 

(-0.45 ± 1.76) x 10" 1 


best fit values for b/a and c/a with errors for CP-even 
higher derivative contribution are given as 

b/a = (1.04 ± 0.43) x 10" 4 GeV" 2 (41) 

c/a = (0.34 ± 1.07) x 10" 4 GeV" 2 (42) 

The la and 2cr contours for b/a vs c/a for Higgs with CP- 
even higher derivative contribution is shown in Fig. 11. 

The values of observables for Higgs in CP-even-odd 
scenario are tabulated in Table VII. It should be noted 
that value of the observable t\ is large due to non zero 
CP-odd admixture for the CP-even-odd scenario. 

At 3000 fb -1 the best fit values of b/a, c/a, 5 for the 
CP-even-odd scenario are given as 

b/a = (0.59 ± 0.40) x 10" 4 GeV" 2 (43) 

c/a = (2.10 ± 0.87) x 10" 4 GeV -2 (44) 


FIG. 11. c/a vs b/a la (green) and 2a (yellow) contours for 
Higgs with CP-even higher derivative contribution at 3000 
fb -1 . The best fit value of (b/a, c/a) is shown by the block 
dot. The values with which data are generated (b/a = 1.25 x 
10” 4 , c/a = 0) is shown by the V. The cross-hair corresponds 
to b = c = 0. 

TABLE VII. The values of the observables for CP-even-odd 
scenario with respective errors at 3000 fb -1 . 


Observables 

Values with errors 

£2 

-0.13 ±0.04 

ti 

0.02 ± 0.03 

Ul 

0.11 ±0.18 

Ul 

(0.27 ± 1.89) x 10" 1 

Vi 

(0.63 ± 1.82) x 10" 1 

Vi 

(0.83 ±1.78) x 10" 1 


5 = 0.57 ± 1.33 in radian. (45) 

The ltr and 2er contours for b/a vs c/a and 8 vs c/a 
are shown in Fig. 12 and Fig. 13 respectively. 

C. Momentum dependence of Higgs couplings 

The vertex factors a, b and c written in Eq. (4) are 
in general momentum dependent but in SM they have 
no momentum dependence. Thus it is essential to ver¬ 
ify their momentum independence to establish H as the 
SM Higgs. To achieve this, one has to measure the mo¬ 
mentum dependence of a, b and c in different momentum 
regions. In LHC one only measures the ratios of cou¬ 
plings i.e. b/a and c/a. However one can measure the 
values of b/a and c/a in different g|(rn 34 ) regions. This 
will allow us to check the momentum dependence of b/a 
and c/a and finally a at 14 TeV 300 fb -1 LHC run. 

Integrating Eq. (11), (12), (13) over qf we get unian- 
gular distribution in n-th bin as follows 

1 dr _ 1 

F n d cos 0 ± 2 


Ti cosdi ± T// -P 2 (cos0i), (46) 
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TABLE VIII. Three y/q. J = m34 bins and corresponding num¬ 
ber of events in each bins for momentum dependence measure¬ 
ments at 14 TeV 300 fb -1 LHC run. 


Bin No. 

Bin range \/gf = 77734 

Number of events 

Bin 1 

12.00 GeV < m 3 4 < 29.00 GeV 

210 

Bin 2 

29.00 GeV < m.34 < 46.00 GeV 

187 

Bin 3 

46.00 GeV < m.34 < 80.00 GeV 

52 


tex factor ‘a’ in these 3 bins, we rewrite the decay width 
in n-th bin r n as follows 

T n = a 2 n T' n (b/a,c/a). (49) 


FIG. 12. c/a vs b/a la (green) and 2<r (yellow) contours for 
CP-even-odd scenario Higgs at 3000 fb _1 . The best fit value 
of (b/a, c/a) is shown by the block dot. The values with which 
data are generated ( b/a = 1.06 x 10 -4 , c/a = 1.89 x 10 -4 ) is 
shown by the ‘> 1 =’. The cross-hair corresponds to b = c = 0. 



c/a (10 -4 GeV -2 ) 

FIG. 13. 5 vs c/a la (green) and 2cr (yellow) contours 

for CP-even-odd scenario Higgs at 3000 fb _1 . The best fit 
values (b/a,c/a) is shown by the block dot. The cross-hair 
corresponds to 8 = c = 0. The value with which data are 
generated is shown by the V. 


where a n is the value of a in n-th bin. T' n is obtained by 
dividing r„ by a 2 and making it a function of b/a and 
c/a only. We can calculate the values of V n and its errors 
in different bins by substituting the values of b/a and c/a 
from Table IX. 


TABLE IX. b/a and c/a with corresponding errors in 3 bins 
in GeV 2 


Bin No. 

b/a in 10" 4 GeV- 2 

c/a in 10 4 GeV 2 

Bin 1 

Bin 2 

Bin 3 

0.02 ± 0.44 

0.39 ± 0.79 

1.35 ±2.09 

0.30 ±0.57 

0.60 ± 0.89 

1.56 ±2.24 


Fit values of b/a and c/a have relatively larger errors 
in last bin due to small number of events and both b/a 
and c/a are consistent with zero in each 77134 bins. 

For resonant production of Higgs we can factor out 
the production cross-section and r n is proportional to 
the number of events N n in each of the three 77734 bins. 
If Tj, and W; , Nj are the value of decay widths and 
number of events in i and j- th bin respectively, then 


Ti 


Nj 


(50) 


< 47 > 

— = — + UV/ cos c/y + U// cos 2 <j> 

T n dcjy 2 it 1 2 

+ V" sin (j> + V 2 sin 2<jy, (48) 

where 77\ 77' , U/', IA% VJ and V” are observables in 77 -th 
bin integrated over range in that bin. 

In SM the vertex factors are momentum independent 
and have constant values. This implies that a must have 
a constant value 1 and b, c should be zero in each 77734 
bin. As a result the values b/a , c/a will be zero in each 
77734 bin. I 11 this subsection we bin the events and find 
the asymmetries T / L , 77\ U /‘, U£ and utilize them to find 
out b/a and c/a in three different bins. 

Now we will study the momentum dependence of ver- 


will hold between any two bins. Furthermore putting 
and r ? expressions as written in Eq. (49) into Eq.(50) 
one finds the following relationship 


Oi 



(51) 


If a is independent of momentum the ratio will al¬ 
ways be 1 for any two bins. We tabulate all possible 
in Table. X for this three bins with their corresponding 
errors. One can similarly perform the same analysis at 
3000 fb _1 . At 3000 fb _1 the errors in b/a and c/a will 
be reduced and with the enhanced statistics, one may in 
principle have enough events to increase the number of 
bins to check momentum dependence of Higgs couplings. 
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TABLE X. Ratio i\j between different bins. 


Ratios (rij) 

values with errors 

ri2 

1.06 ±0.28 

r 13 

0.82 ±0.38 

V23 

0.75 ±0.39 


These Ratios are consistent with SM (i.e ratio = 1) 
as they should. The ratio 7-13 and have larger errors 
due to small statistics. These results will be improved by 
the bin size optimization and with larger statistics. 

IV. COMPARISON OF PRECISION 
MEASUREMENT BETWEEN 33 TEV LHC AND 

ILC 

After HL-LHC (High Luminosity LHC), LHC energy 
could be upgraded to run at a center of mass energy of 
1 /s = 33 TeV. In this section we compare the precision 
measurements of H coupling between 33 TeV LHC and 
ILC. We have followed the same cut based analysis as 
14 TeV 3000 fb~ Y machine. At 33 TeV 3000 fb " 1 the 
number of events are given in the cut flow Table XI. The 


Cuts 

Number of events 

Selection cuts 

6523 

50 GeV < mi 2 < 106 GeV 

6362 

12 GeV < 77134 <115 GeV 

5935 

115 GeV < mu < 130 GeV 

5852 


TABLE XI. Effect of the sequential cuts on the simulated 
Signal for 33 TeV 3000 fb -1 LHC upgrade. The fc-factor for 
signal is 1.95. 



b/a (10 -4 GeV -2 ) 

FIG. 14. c/a vs b/a la (green) and 2a (yellow) contours for 
SM Higgs at 33 TeV 3000 fb -1 LHC run. The best fit values 
(b/a, c/a ) is shown by the block dot. The V corresponds to 
b = c = 0. 

can find out the partial width T(H —> ZZ) 

T(H^ZZ) 

total BR(H -A ZZ) 1 ] 

For example the expected precision for CP-odd anoma¬ 
lous coupling c at ILC [64] is 7 x 10 -4 to 8 x 10 -6 which 
is roughly around the loop induced CP-odd contribution. 
We have shown that precision in the measurement of c/a 
can be 6 x 10 ~ 5 for 33 TeV LHC using angular asymme¬ 
tries. However at ILC the measurement of Higgs decay 
width will allow us to extract the absolute values a, b, c 
which is beyond the scope of LHC. 


V. CONCLUSION 


best fit values of b/a and c/a for SM Higgs are 

b/a = (2.64 ±1.60) x 1(T 5 GeV -2 (52) 

c/a = (1.07 ± 6.25) x 10 ~ 5 GeV -2 (53) 

In our first attempt, we find that a precision of 10 -5 
is achievable at the 33 TeV 3000 fb -1 LHC, for CP-odd 
admixture and CP-even higher derivative contribution. 
In future, a more detailed analysis may improve the result 
and decrease of errors in b/a and c/a. Although the 
precision achieved is high, LHC will only measure the 
ratios of Higgs couplings even at 33 TeV. 

ILC will provide an independent test in measuring 
HZZ couplings subject to its proposed implementation 
and will also offer invaluable probe to such couplings. 
At LHC one measures a ■ BR(H -A ZZ) however ILC 
will measure the branching ratio BR(H -A ZZ) by mea¬ 
suring the inclusive cross section azH for the process 
e + e _ —> ZH. This inclusive measurement of cross sec¬ 
tion alone will probe HZZ couplings to 1.3%[63]. Iden¬ 
tifying a Z boson in recoil against the Higgs boson one 


We demonstrate that angular asymmetries will pro¬ 
vide a strong and efficient tool to probe Higgs couplings 
in high luminosity future LHC runs. With the increased 
statistics at 14 TeV run, LHC will enter into precision era 
and angular analysis will offer a step by step methodol¬ 
ogy to study the Higgs couplings. Angular asymmetries 
can be utilized to probe the Higgs couplings and to dis¬ 
entangle its exact CP property in the next LHC run. We 
benchmark our observables for SM, CP-odd admixture, 
CP-even higher derivative contribution and finally CP- 
odd admixture with Higher derivative CP-even coupling. 
We perform the analysis for two different luminosities at 
300 fb -1 as well as 3000 fb -1 and study the precision 
with which angular analysis probe HZZ couplings. The 
study of the momentum dependence of the Higgs cou¬ 
plings would be a significant step forward in establish¬ 
ing its SM nature, since in the SM Higgs couplings do 
not have any momentum dependence. At 14 TeV LHC 
run with the improved statistics, we present how one can 
examine the momentum dependence of the Higgs cou¬ 
plings in different momentum regions. We have further 
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discussed what precision LHC can achieve in the mea¬ 
surement of the Higgs couplings to Z boson by angular 
analysis at 33 TeV. A comparison has also been made be¬ 
tween precision measurement of 33 TeV LHC and ILC. 
Angular analysis will be a powerful technique to decipher 
the CP properties of Higgs couplings at 14 TeV LHC run 
and will open up a new domain of precision measurement. 
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Appendix A 


|Ai — A 2 | ^ J and the helicity amplitudes are related 
as Aa 2 ,Ai = (— 1 ) j A_Ai,-a 2 - AVe have three orthogonal 

helicity amplitudes .Too, -4++ and A _However one 

can write three helicity amplitudes in transversity basis 
as: 


Aoo 

(A2) 

-j=(A ++ + A —) 

(A3) 

-±=(A ++ + A-). 

(A4) 


Helicity fractions Fl, F\\ and F±_ are defined as 

F a = , where A £ {L, ||, _L}. 

\/|- 4 l | 2 + |. 4 || | + 1-4 j _| 2 

Moreover n = with v e = 2 1 3 e — 4e t sin 2 9w and 

V £ 

a e = 2 / 3 ^. 

Tf is defined as 


The amplitudes for the process H with spin J (spin 0) 
decays to two Z(spin 1) boson with spin projections 
along z axis Ai and A 2 is [67, 68 ] 

/ p T + 1 \ 3 

Ai.Aa) = (■ 

(Al) 

where @j\ is the Wigner-D function, A = | Ai — A 2 1 and 
-4 ai,a 2 are the helicity amplitudes with Ai^ £ {± 1 , 0 }, 
J = |J|. Angular momentum conservation implies |A| = 

I 


r f = 


dr 

dql 


Af ^|.4l | 2 + |-4|| | + |-4j_| 2 ^ 5 (A5) 


Br i 


X 


with M cos 2 g w jgr- Mz ^ q 2_ M 2 z y +M 2 T 2 z ^ ■ 1 z 1S 

the decay width of Z, Br« is the branching fraction for 
the decay of Z to two massless leptons. However we as¬ 
sumed narrow width approximation for on-shcll Z\ boson 
in Sec.II. In Sec.Ill we have implemented the the cut flow 
table while integrating over q\ and to find the expres¬ 
sions for 7{, T 2 , U\ , U 2 , Vi and V 2 . 

The expressions for 71, 72 , U\, U 2l Vi and V 2 are 


Ti 

r 2 

u x 

u 2 

Vi 

v 2 


1.32 x 10 ~ 9 y 

5.57 x 10- 8 ± 2.61 x 10~ 8 :r + 3.98 x 10" 9 x 2 ± 1.60 x lO^y 2 
-9.65 x 10 -9 + 4.00 x 10 _lo x 2 + 4.00 x lO” 10 ^ 2 

5.57 x 10~ 8 + 2.61 x 10~ 8 x + 3.98 x 10“ 9 x 2 ± 1.60 x 10- 10 y 2 
-1.20 x 10 -9 - 6.33 x 10~ 10 a; ± 7.11 x 10" n x 2 

5.57 x 10- 8 ± 2.61 x 10~ 8 x + 3.98 x 10~ 9 x 2 ± 1.60 x 10- 10 y 2 
6.06 x 10~ 9 + 4.35 x 10" 9 x ± 7.97 x 10" lo x 2 ± 4.00 x 10 ~ 10 y 2 

5.57 x 10- 8 + 2.61 x 10- 8 z + 3.98 x 10~ 9 a; 2 + 1.60 x 10- W y 2 
3.11 x 10 -lo y 

5.57 x 10 - 8 ± 2.61 x 10~ 8 x + 3.98 x 10~ 9 x 2 ± 1.60 x lO^y 2 
-2.92 x 10~ 9 y 

5.57 x 10 ~ 8 + 2.61 x 10~ 8 x + 3.98 x 10“ 9 x 2 ± 1.60 x 10- 10 y 2 


where x = % x (lOOGev ) 2 and y = ^ x (lOOGev ) 2 


(A6) 

(A7) 

(AS) 

(A9) 

(A10) 

(AH) 
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